Quantum turbulent velocity statistics and quasiclassical limit 
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Two research groups have measured turbulent velocity statistics in superfluid helium using dif- 
ferent techniques. The results were in conflict: one experiment revealed Gaussian distributions 
(as observed in ordinary turbulence), the other experiment determined power-laws. To solve the 
apparent puzzle, we numerically model quantum turbulence as a tangle of vortex filaments, and 
conclude that there is no contradiction between the two experiments. The transition from Gaussian 
to power-law arises from the different length scales which are probed using the two techniques. We 
find that the average distance between the quantum vortices marks the separation between quantum 
and quasi-classical length scales. 



Turbulence in superfluid helium (''He) [H Q consists 
of a tangle of discrete, thin vortex filaments, each car- 
rying one quantum of circulation k = h/m (where h is 
Planck's constant and m the mass of one atom). This 
'quantum turbulence' 0, 0] is also studied in superfluid 
■^He 0,0] and in atomic Bose-Einstein condensates Q- 
The absence of viscosity and the discrete nature of the 
vorticity distinguish quantum turbulence from ordinary 
(classical) turbulence. Nevertheless, experiments have 
also shown remarkable similarities: the same pressure 
drop along and pipes 0] , drag on a moving body Q , vor- 
ticity decay [lo| , and distribution of kinetic energy over 



the length scales (Kolmogorov energy spectrum) [11|, [12 1 . 
The relation between quantum turbulence and classical 
turbulence has thus become the focus of attention. 

This report is concerned with the statistical proper- 
ties of the components of the turbulent superfluid ve- 
locity field V. In ordinary turbulence, experiments [isj 
and numerical simulations [H, 15| confirm that velocity 
statistics are Gaussian. The usual interpretation of this 
property is the following [16[. At any time t, the velocity 
field V at the point r is determined by the vorticity field 
u; = V X V via the Biot-Savart law [17| 
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Thus, if the point r is surrounded by many, randomly 
oriented vortical structures, Gaussianity results from the 
application of the Central Limit Theorem. 

In the last few years, the development of particle 
tracking visualization techniques suitable for liquid he- 
lium [isl [Toj has opened the way to the measurement 
of quantum turbulence statistics. Using solid hydrogen 
tracers, Paoletti et al. at the University of Maryland 
deduced that, in superfluid helium, the velocity compo- 
nents of the turbulent velocity follow power-law distri- 
butions. Large velocity fluctuations are thus relatively 
more frequent than in ordinary turbulence, where these 
distributions are Gaussian [2l|. The authors suggested 
that power-laws distributions arise from vortex reconnec- 
tions, clearly high-velocity events. 



The experiment of Paoletti et al. motivated theoretical 
work on the problem. White et al. [22| solved numerically 
the Gross-Pitaevskii equation for a Bose-Einstein con- 
densate and determined velocity distributions in a vari- 
ety of quantum vortex systems at zero temperature: two- 
dimensional and three-dimensional trapped atomic con- 
densates, three-dimensional homogeneous condensates, 
Onsager's two-dimensional vortex gas. In all systems, 
they found power-law distributions, in agreement with 
Paoletti et al. [l^, even in the absence of vortex recon- 
nections. Using results from previous work by Min and 
Leonard uSl on the numerical analysis of vortex methods 
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(later confirmed by Weiss et al. |2J]), White et al. 
concluded that the non-Gaussianity arises from the sin- 
gular nature of the quantum vorticity, which causes an 
anomalously slow convergence to Gaussian behaviour. 

Further theoretical work at non-zero temperatures 
seemed to strengthen this conclusion. Baggaley and 
Barenghi j25| used the vortex filament model to numer- 
ically study intense vortex tangles. They confirmed the 
Kolmogorov spectrum and the power-laws distributions 
of the velocity components. The vortex filament model 
was also used by Adachi and Tsubota [IB]: who found 
power-law distributions in counterflow turbulence (quan- 
tum turbulence driven by a heat flow) 

Salort et 



In a recent experiment, Salort et al. [271 at CNRS 
Grenoble measured the turbulent velocity in a superfluid 
wind tunnel using a small Pitot tube. They confirmed 
the Kolmogorov spectrum, but reported Gaussian veloc- 
ity distributions, in apparent contradiction with Paoletti 
et al. [iO] and with the numerical simulations 2h.^ 2Gi, ■ 
To solve this puzzle, firstly we consider the size of the 
probe relative to the distance between vortices in the two 
experiments. The radius a of the tracer particles used 
by Paoletti et al. was a w 10~^ cm. The turbulence 
was created by a heat flux Q which drove a counterflow 
current (28| Vns — l^n — ^sl = Q I {PsST) between super- 
fluid and normal fluid velocity components w„ and Vs, 
where ps is the superfluid density and S the specific en- 
tropy. For example, at T = 1.9 K the largest heat fiux 
was Q = 0.17 W/cm^, corresponding to Vns = 1.5 cm/s, 
with p, = 0.0831 g/cm3 and S = 0.709 J/(gK). Using 
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the relation L = 7^w,^js with 7 « 140 s/cm^ as reported 
in Ref. [2^, one finds £ w 4.7 x 10"'^ cm, where L is 
the vortex fine density (vortex length per unit volume) 
and (. ~ Paoletti et al. [20[ took their velocity 

measurements during the decay of the turbulence, after 
switching off the heater, so the actual inter- vortex dis- 
tance must have been larger than the above value of £. 
We conclude that, in the experiment of Paoletti et al. (20l |. 
a/£ ^ 1; this means that, although the tracer particles 
must have altered in some measure the motion of the 
vortex lines, the velocity field was probed at very small 
length scales, between £ and a. 



In the experiment of Salort et al. [27|, the nozzle of 
the Pitot tube had diameter 0.06 cm. The authors esti- 
mated [2^ that the inter-vortex separation was between 
£ sa 5 X lO"'^ and 20 x lO'^cm. We conclude that in 
this experiment a/i ^ 1 (even without considering that 
the effective resolution in the streamwise direction, set 
by time dynamics of the probe, was perhaps much larger 
than a|29|). In summary, the experiment of Paoletti 
et al. [20| probed the velocity field between vortices, at 
length scales less than £, whereas the experiment of Sa- 
lort et al. [i^] detected motion at scales larger than 
containing many vortex lines. 

Secondly, we perform the following numerical calcu- 
lation. In both experiments the length scales involved 
were much larger than the vortex core radius, oq = 
10~^ cm, so it is appropriate to use the vortex filament 
model [30[ • We assume that the vorticity is concentrated 
in space curves s(^, i) around which the circulation is 
K = lO^^cm^/s. Eq. ([T]) reduces to 
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where the line integral extends over all vortex lines. The 
numerical techniques to discretize the vortex lines, to 
compute the time evolution, to de-singularize the Biot- 
Savart integral and to algorithmically perform reconnec- 
tions when two vortex lines become close to each other 
are described in the literature [2^, [s^ or in our previous 
papers pH 31 1. The computational domain is a peri- 
odic box of size D = 0.075 cm. The distance between 
discretization points along the vortex filaments is held 
at approximately 6/2 where 5 ~ 0.001 cm. To speed- 
up the evaluation of Biot-Savart integrals we use a tree 
method [s^] with critical opening angle 9 = 0.4. For 
the sake of simplicity, wc assume that the temperature 
is T = 0. This means that no external forcing is needed 
to sustain the turbulence, because the total kinetic en- 
ergy is conserved during the evolution (at least for the 
time scale of interest here; the finite discretization and 
the reconnection algorithm introduce small energy losses 
which model phonon emission as described in Ref. (3l|). 

The initial condition at time t = Q consists of 100 
straight, randomly oriented, vortex filaments discretizcd 




FIG. 1. Snapshot of the turbulent vortex tangle at t = 0.1 s 
after the vortex line density has saturated, the box size D = 
0.075 cm. 



over TV = 10^ Lagrangian points Sj (j = 1, • • • N). During 
the evolution the vortex lines interact, become curved, re- 
connect, and a vortex tangle is quickly formed. We find 
that the vortex line density L initially grows, then satu- 
rates to the value L sa 2.04 x 10'* cm"^. We check that 
the average curvature saturates too. Fig. [T] shows the 
vortex tangle when we stop the calculation at t = 0.1 s. 
We examine the energy spectrum E{k), defined by 
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E{k)dk, 
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where E is the total energy per unit mass. Fig. [2] 
shows that most of the energy is contained in the large 
scales, and that the spectrum is qualitatively consistent 
with the Kolmogorov scaling E^^ ^ k^^/'^ in the region 
fc <JC 27r/^ « 900 cm~^, in agreement with experiments 
[ill [1^ and with previous numerical calculations [25 , 33 - 
[35[ . The shallow spectrum at large values of k is also con- 
sistent with existing similar calculations [s^]- Secondly, 
we examine the velocity components and compute their 
probability density functions (normalized histograms, or 
PDF for short). Rather than sampling the values at 
points of the domain, we compute averages over regions 
of size A, varying A. This is done by calculating the 
velocity field, induced by the vortices, on a 128'^ Carte- 
sian mesh, and then averaging cells over the appropriate 
length scale. The result is shown in Fig. [31 where we 
compare the data against their Gaussian fits 
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FIG. 2. (Color online) Energy spectrum Ek (arbitrary units) 
vs wavenumber k (cm~^) corresponding to Fig.[T] The dashed 
line is the Kolmogorov scaling Ek ~ k~^^^; the solid grey line 
shows wavenumber corresponding to the intervortex spacing 
ke f» 900 cm"\ 



where /i w and a is the standard deviation. It is 
apparent that, if A < £, the PDFs are not Gaussian; 
by fitting PDF(ui) - i;^'' we obtain b = 3.2, 3.2 and 
3.1 for i = x,y,z respectively, in agreement with the 
power-law statistics found experimentally by Paoletti et 
al. [2^ and numerically by White et al [22j, Baggaley 
and Barenghi [2^, and Adachi and Tsubota [2Q|. If 
A > £, however, the distributions are Gaussian, and we 
recover the same (classical) statistics measured by Salort 



et al. [27 1 



In conclusion, the different velocity statistics observed 
by the Maryland and Grenoble experiments are consis- 
tent with each other. Taken together, these experiment 
and this work support the interpretation that, at scales 
larger than £, quantum turbulence exhibits quasi-classical 
behaviour (in terms of both energy spectrum and veloc- 
ity statistics), whereas at scales smaller than £ (but still 
orders of magnitude larger than the quantum coherence 
length ao), the discrete nature of quantized vorticity af- 
fects the distribution of energy over the length scales and 
the frequency of high- velocity events. 
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by the Leverhulme Trust (Grants F/00125/AH and RPG- 
097). 




FIG. 3. (Golor online) Probability density functions (PDF) 
of turbulent velocity components Vi (i = x, y, z) vs Vi/ai cal- 
culated by averaging over regions of size A = 2^ (top), l\ = I 
(middle) and A = l/Q (bottom). (Green) asterisks, (blue) 
circles and (red) triangles refer respectively to i = x, i — y 
and i = z components. The solid lines ares the Gaussian fits, 
defined by Eq. H where a = 0.041 (top), 0.048 (middle) and 
0.082 cm/s (bottom) respectively. 
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